Black Holes in Higher Dimensional Gravity Theory with Quadratic in Curvature 

Corrections 



Valeri P. Frolov* and Ilya L. Shapiro^ 

* Theoretical Physics Institute, 
Department of Physics, University of Alberta 
Edmonton, AB, Canada, T6G 2J1 
E-mail: frolov@phys.ualberta.ca and 

1 Departamento de Fisica, ICE, 
Universidade Federal de Juiz de Fora 
CEP: 36036-330, Juiz de Fora, MG, Brazil 
E-mail: shapiro@fisica.ufjf.br 
(Dated: August 12, 2009) 

Static spherically symmetric black holes are discussed in the framework of higher dimensional 
gravity with quadratic in curvature terms. Such terms naturally arise as a result of quantum 
corrections induced by quantum fields propagating in the gravitational background. We focus our 
attention on the correction of the form C 2 = C a (3~ l &C otli ~ 1& . The Gauss-Bonnet equation in four- 
dimensional (4D) spacetime enables one to reduce this term in the action to the terms quadratic in 
the Ricci tensor and scalar curvature. As a result the Schwarzschild solution which is Ricci flat will 
be also a solution of the theory with the Weyl scalar C 2 correction. An important new feature of the 
spaces with dimension D > 4 is that in the presence of the Weyl curvature-squared term a solution 
necessary differs from the corresponding 'classical' vacuum Tangherlini metric. This difference is 
related to the presence of secondary or induced hair. We explore how the Tangherlini solution is 
modified by 'quantum corrections', assuming that the gravitational radius ro is much larger than 
the scale of the quantum corrections. We also demonstrated that finding a general solution beyond 
the perturbation method can be reduced to solving a single third order ODE (master equation). 

PACS numbers: 04.50.-h, 04.60.Bc, 04.70.-s, 97.60.Lf Alberta.-Thy-08-09 



I. INTRODUCTION 

The Einstein-Hilbert action is the simplest possible 
gravitational action. It can be easily generalized to the 
higher dimensional models of gravity 

W=^Jd D x^gR. (l) 

Black hole solutions with the spherical topology of the 
horizon in this theory are well known now. The most 
general solution (in the presence of a cosmological con- 
stant) is a so- called Kerr-NUT-(A)dS spacetime [l|. In 
the absence of rotation and for vanishing NUT parame- 
ters and the cosmological constant, this solution reduces 
to the Tangherlini metrics 

ds 2 = -Adt 2 + ^ + r 2 duo 2 n , A =l-(^y~\ (2) 

For a Z?-dimensional spacetime, n = D — 2, and dtu 2 is 
the metric of a unit round sphere S n . 

Modifications of the Einstein-Hilbert action are com- 
monly considered in the modern literature. Besides the 
needs of a purely phcnomcnological description of grav- 
ity in different models, there exists a more fundamen- 
tal reason for considering a generalized Einstein gravity. 



It is well known that quantum corrections in a gravita- 
tional field can be described by using the DeWitt's ef- 
fective action formalism (3[. In a general case, such an 
effective action contains higher in curvature corrections 
to the Einstein-Hilbert Lagrangian, as well as non-local 
contributions (see, e.g., review [^| and references therein). 
Higher in curvature corrections were used for example in 
the Starobinsky model of inflation The higher cur- 
vature terms naturally arise in the string theory (see e.g. 
@ and references therein). One of the most interest- 
ing aspects connected to the introduction of the higher 
in curvature effective actions is a possibility to resolve 
a curvature singularity problem. Such singularities, ac- 
cording to the Penrose and Hawking 0, H @ theorems 
always exist in cosmology (initial singularity) and inside 
black holes (final singularity). It is well known that the 
interior of (non-rotating) black hole in the vicinity of the 
singularity is similar to a contracting anisotropic universe 
and it exhibits the Kasner type behavior with infinitely 
growing curvature invariants. In the models where the 
curvature is limited by some (say Planckian) values, one 
may expect newly born universes instead of the singu- 
larity formation (see, e.g. [HI EH)- For more general 
discussion of the quantum effects in black holes see e.g. 
the book [l^] and references therein. 

Let us notice that the consistent consideration of quan- 



2 



turn effects in gravity requires an introduction of the 
higher derivative terms, even at the semiclassical level 

(see m m for the introduction and 0] for the recent 

review of semiclassical approach) . One of the first papers 
where the quantum curvature corrections were considered 
in the gravitational collapse model was [15j |. In particu- 
lar, it was shown that a collapse of a null shell with the 
mass smaller that the Planckian mass does not create a 
black hole. In other words, in this problem there exists 
a mass gap for the black hole formation. In the modern 
language it means that the black hole formation in the 
gravitational collapse of null shells is a first order phase 
transition. Later, black holes in the theories with higher 
in curvature terms were considered, e. g. in [l6l |. 

In the four-dimensional (41?) case the lowest in cur- 
vature quadratic corrections to the Einstein gravity do 
not modify the Schwarzschild solution. The reason is 
the following. In AD case there exists the Gauss-Bonnet 
relation 

C abcd C abcd -2^R ab R ab -^R 2 ^j = topological term. (3) 

As a consequence, the general AD action with quadratic 
in curvature corrections can be always written as 

W = -!- / d 4 x,/=g~ (R + aR ab R ab + bR 2 ) . (4) 

167T J V ' 

It is clear that any vacuum solution of the Einstein equa- 
tions, Rob = is, at the same time, a solution of the the- 
ory ((4]). For this reason, in order to study the quantum 
gravity corrections we need to use higher in curvature 
Lagrangians. One of the interesting option is to consider 
AD theories with the corrections of the form f(C 2 ), which 
in many respects are similar to the f(R) theories in the 
cosmological models. 

In this paper we use another approach. We consider 
higher dimensional spherically symmetric black holes 
with quadratic in the curvature corrections. Since the 
Gauss-Bonnet argument does nor work directly in the 
higher dimensions, the quadratic in the Weyl tensor cor- 
rections necessarily modify the background Tangherlini 
solution. This results in the creation of what is called 
secondary or induced hair fl7l . [l8| . In the present pa- 
per we focus our attention on the study of these 'hair'. 
Namely, we consider static spherically symmetric vacuum 
solutions of the gravitational theory with C 2 corrections. 
Using the reduced action approach we derive equations 
for such spacetimes (sections II and III) . In section IV we 
demonstrate that in the special (seminull) coordinated 
these equations can be reduced to a single third order 
ordinary differential equation (ODE) for a quantity, con- 
nected to the Weyl curvature. Linearized version of the 
field equations and their static solutions are discussed in 
section V. In order to obtain corrections to the Tangher- 
lini metric describing the induced hair generated by the 



C 2 correction, we use the iteration procedure developed 
in section VI. This section also contains the numerical 
results. Specific properties of the AD case are discussed 
in section VII. Section VIII contains discussion and lists 
some of the perspective problems. 

Throughout the paper we use the sign conventions 
adopted in [lj5[ . 



II. REDUCED ACTION 

In a D-dimensional spacetime with metric 

ds 2 = g^dx^dx 1 ' (5) 
we consider a gravitational action of the form 

w - il/^^-fO. < 6 > 

where C 2 = C a0lS C a01& . (7) 

Here R and C a f3ys are the Ricci scalar and the Weyl ten- 
sor, respectively. The constant a has the dimensionality 
of [length] 2 , and, for physically realistic case, it is posi- 
tive. The reason is that this sign corresponds to the pos- 
itively defined theory, e.g., to the positive energy of the 
massless mode (graviton) in the flat space limit. This re- 
quirement is also relevant because it provides the correct 
low-energy limit of the theory [3, |2Q| • The constant a 
can be interpreted as a square of the characteristic funda- 
mental scale I for an appropriate quantum gravity theory. 
In a 'standard quantum gravity' I is of the order of the 
Planckian length. In the popular recently models with 
large extra dimensions the corresponding to I fundamen- 
tal energy scale is of order of TeV. We shall study spher- 
ically symmetric black hole solutions. The characteristic 
scale of such a solution is its gravitational radius r . This 
dimensional parameter can be chosen as a general metric 
scale. After this in the theory © there remains only one 
essential dimensionless parameter 1/rQ (or, equivalently, 
a/r 2 ,). One can expect that when ro ^> I, the Weyl scalar 
(Weyl curvature-squared) term plays the role of a small 
correction. We shall focus on this case later on, but at 
the initial part of the work there is no need to make ad- 
ditional assumptions about this parameter. 

Our purpose is to explore static spherically symmet- 
ric solutions of the theory formulated above. The corre- 
sponding metric can be written in the form 

ds 2 = r 2 ds 2 = r 2 (dj 2 + duj 2 n ), (8) 

where n = D — 2 and g?o-> 2 is the metric on a unit sphere 
S n 

dujl = u} ij d0 i d0 j , . . . = 2, . . . , n + 1 . (9) 
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Quantity r is a scalar function on a 2D manifold with 
the metric 



d^ 2 = j a bdx a dx b , a, 6, . 



0,1 



(10) 



For this class of the metrics the action ([6]) can be reduced 
to the 2D action. Because of the spherical symmetry 
neither R nor C 2 depend on the angle variables. Using 
the relation 



1 =r n+ \ 



"7v w 



and integrating over the unit sphere one obtains 



W 

c 



-W 

,™+2 , R 



w 

a „• 



d 2 xy 



-it; 



(11) 
(12) 

(13) 



Here 



27r(" +1 )/ 2 



is a surface area of the n-dimensional unit sphere S n . 

We denote by symbols with bars the objects con- 
structed for the metric ds 2 . Since this metric is a direct 
sum of the two metrics, d'y 2 and du; 2 , the curvature for 
it is also a direct sum of two curvatures and it has the 
following non-vanishing components 



Rabcd = 2 R {lacjbd — Jadjbc) 



Rijkl 

We also have 



(14) 



1 



Rac = 2 R la 



R 



Rik = (n- l)uj lk 
R + n(n-l). (15) 



We denote by symbols with hats the objects con- 
structed for the 2D metric dj 2 . In particular, R and 
C 2 denote the Ricci scalar and the Weyl-square invari- 
ant, calculated for the metric y a b. It is easy to show that 
(see, e.g., (U) 



R = r" 



□r 

R + n(n + l)-2(n+l) — 

r 



(Vr) 2 

- (" + 1 )(«- 2 )t 



(16) 



The conformal invariance of the Weyl tensor implies 
C 2 = r~ A C 2 . (17) 
To calculate C 2 we use a relation 



C 2 =R, 



R 



-R a pR a ^ 



n(n + 1) 



R z 



(18) 



By using (fT4"|) - (fT5")) one can obtain 



Raf3~fS 



R 2 + 2n(n- 1) , 



R: 



l 



a/3 



R 2 + n{n-lf 



R z = 



R + n(n - 1) 



(19) 
(20) 

(21) 



Then we find 



C 2 = 



n- 1 
n + 1 



(R + 2) 2 



Thus the dimensionally reduced action on the 2D man- 
ifold with the metric j a b is 

W = W (c) + W (q) = J d 2 x^j [£ {c) + C (q) ] , (22) 



£ (c) = r n R + n(n - l)r" + n(n + l)r n " 2 (Vr) 2 , (23) 

^ ) =- a( r„ 1 , ) :r 2 (^+2) 2 . m 



2(n + l) 



III. FIELD EQUATIONS 

The reduced action W is a functional of the 2£> metric 
"f a b and of a scalar function r. Its variation has the form 



6W 



d 2 x^l [V8r + g ab 5 lab ] 



(25) 



Consider a coordinate transformation generated by a vec- 
tor field £° 

x a ^x a + C- (26) 
Under this transformation r and 7 a & transform as follows 

S(r = -er, a , Saab = -2V ( „&) . (27) 

Since W is invariant under this transformation one has 
5(:W = 0. Using and (J27J) one obtains the following 
Noether identity 



2V b G b a - VV a r = . 



(28) 



To obtain the field equations for the reduced action we 
use the following relations 



SV~7 = 7;V^fi ab hab ■ 



1 

2 

-(6R)S = Jd 2 V^S Jab 

-j ab as-R ab s . 



v a v b s 



(29) 
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The variation of (12"2")) with respect to r gives 

v = v {c) +v {q) , 

V [c) = nrK -n(n+l)r n - 3 [(n-2)(Vr) 2 
+ 2rUr - (n - 2)r 2 ] , 



V 



(<-/) 



o(n- l)(n- 2) r2 
2{n+l)r n ~ 1 



(30) 

(31) 
(32) 



Here and later we use subscripts c and q for quantities 
connected with 'classical' and 'quantum' parts of the ac- 
tion, respectively. Similarly, the variation of (|22p with 
respect to j a b gives 



/-tab /->ab , r->ab 

y = y {c) + y {q) , 



(33) 



where 



fab 



'(c) 



= nr n - 2 {rV a V b r - 2V a rV h r (34) 
+ \l ab [i n ~ l > 2 - (n~ 3)(Vr) 2 - 2rfir]| 



IV. EXPLICIT FORM OF THE FIELD 
EQUATIONS 

In the absence of the Weyl term in the action one can 
show that any spherically symmetric vacuum solution of 
the higher dimensional Einstein equations is static (a gen- 
eralized Birkhoff 's theorem) . In the presence of the Weyl 
term this is not true anymore. In the present work we 
restrict ourselves by considering static solutions of the 
theory (0), that is a special subclass of its solutions. 

Let us describe the form of the action which is well 
adapted to our problem. We use the coordinate freedom 
r — > r(r) in order to put g rr = 0, and write the metric 
j a b in the form 



A 

d<y 2 = -—e 2C dv 2 + — e c dvdr , 



2 



(40) 



where A = A{r) and C — C(r) are some unknown func- 
tions. The corresponding 'physical' metric is 



ds 2 



-Ae 2C dv 2 



2e c dvdr 



r 2 dtul 



(41) 



and 



r>ab 



a(n — 1) 



(n + 1) 



2-n ts2 



V a V b K + j a \K - UK 



K 2 ) 



(35) 



In the above relations we use, instead of R, the related 
quantity, 



K = r"- 2 (i? + 2), 



(36) 



which is directly connected to the Weyl scalar in the 
'physical' spacetime 



\K\ 



n+1 



r"(C 2 ) 1/2 ■ 



n - 1 

Furthermore, it is easy to show that 



o(n - 2) (n-1) 2 ^ 
4(n+l)r™" 1 a ' 



(37) 



(38) 



This relation allows one to check the validity of the 
Noether identity (|2"8|) for the 'quantum' part of the ac- 
tion. This serves as a good test of the rather long calcu- 
lations required to derive relations (j30 ]) - (f3"5l) . 

To summarize, the vacuum field equations for the re- 
duced action (l22l) are 



V = 0, Gab=Q- 



(39) 



We assume that the spacetime is asymptotically flat. 
This implies that A(r — oo) = 1. The function C(r) 
takes a finite value C at infinity. This value depends on 
the normalization of the advanced time coordinate v. We 
choose its normalization so that C = 0. For this choice 
£ = d Vl which is a Killing vector, has a canonical normal- 
ization at the infinity: £ 2 = — 1. 

The event horizon, if it exists, is located at the sur- 
face where A = 0. The metric (|4"Tj) in the coordinates 
(v, r) is regular at this surface, provided C is regular 
there. In particular, this property enables one easily to 
consider modifications of the Schwarzschild black hole 
solution due to the quantum corrections, in a case when 
the latter are small and can be treated as perturbations 
(see section 4). In the presence of a black hole, its event 
horizon is located at the point th where the function A 
vanishes. The surface gravity k of a static black hole is 
defined as 



1/2 



(42) 



where £ = d v is the Killing vector. For the metric 
this expression can be rewritten in a more simple form, 



1 



K=^-{e u A')\ 

2 v / | t—Th 



(43) 



The Hawking temperature of such a black hole is Tjj = 
2tt/k. For the unperturbed metric, when a = 0, one has 



K 



n-1 
2r 



(44) 
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Calculations using GRTensor program allows one to 
obtain the following relations. The non- vanishing com- 
ponents of the Christoffel symbols are 



r, 



r, 



v, 



AVe 



= C - -. 

r 



where 



V = T-(A'r + 2AC'r - 2 A) e c . 
2r 



(45) 



(46) 



One also has Q vv = Q vr , due to algebraic properties of the 
metric f|40[) and corresponding components of the Ricci 
tensor. Finally, it proves useful to introduce a parameter 
p, related to a as 



P 



2(n + 1) 



Notice that according to our assumption this parameter 
is negative. The non-vanishing field equations can be 
presented in the form 



(n - l)p Ar 2 AK" + 2r(rA' + 2A)K' - 4K 
+nr n [rA' + (n - 1)(A - 1)] = , 



K 



r n — 2 



(47) 



nr n C - 2(n - l)p[rK" - rC'K' + 2K'] = , (48) 

(n - l)(n - 2) P K 2 r 2i - l ~ n *> - n [(n - 2)(n + 1)(A - 1) 
+2(n + l)r{AC + A') - K r 2 ""] = , (49) 



K 



_n— 2 



[3r 2 A'C - 2 - 2r(A' + AC') 



+2A + r 2 A" + 2Ar 2 {C" + C" )] 



(50) 

The equations (|47|) and ([48]) follow from (vr) and (rr) 
components of the equation Q a b = 0, while (|49|) follows 
from the equation V = 0. The last equation is the defi- 
nition of K . 

A remarkable property of the system of equations (1471) 
(US]) is that they can be reduced to a single third order 
ODE. Wc notice first that this system does not contain 
the function C(r), but only its derivatives. The ([4*5]) 
enables one to express C as a function of K and its 
derivatives 



c = 



2(n- l)p(rK" + 2K') 
nr n + 2{n — l)prK' 



(51) 



After substituting this expression into (|47|) and (|49|) , one 
obtains two linear equations for A and A'. Solving these 
equation one determines A and A 1 as functions of K and 
its derivatives up to the second order. One can present 
the expression for A in the form 



A = 



SiK' + S a 



n{T 2 K" + T 1 K> + T ) ' 



(52) 



where 

Si = -2p(n - 1) x 



i(n - 2)(n + l)r - p(ri 3 - 3n + 2) 



X 2 



8(n 2 - l)pK - nr 2 K 



p(n- l)(ri + 4)is: 2 



T 2 = 4(n 2 -l)pr 2 , 
T = n 2 (n + l)r n . 



Ti = -2p{n 2 - l)(n- 2)r, 



A similar expression is valid for ^4'. Taking a deriva- 
tive of A given by (|52|) and putting it equal to A' , one 
arrives at a single ordinary differential equation of the 
third order for K . We call it master equation. We have 
obtained this equation using the MAPLE program but 
do not reproduce it here since it is too bulky. 



LINEARIZED EQUATIONS AND THEIR 
SOLUTIONS 



Since the equations of our interest are rather compli- 
cated, it is worthwhile to consider first their linearized 
form. This exercise is especially useful, because it helps 
to control the correct classical (p — ► 0) limit in the region 
far from the center r = 0. Let us denote 



A = l + a{r) , C = P(r) , K = v{r) 



(53) 



We regard the functions a, /3, and v as small perturba- 
tions and hence, in the course of our calculations, we shall 
keep only those terms which are linear in these functions. 
The linearization of equations (|4"T)) - (f4"9"|) gives 



n(n — l)a + nrct 
-A{n-l)pr- n (v- 



(54) 



TV - T 2 V"\ 



= 0, 



2(n - l)p(ru" + 2v') - nr n = , (55) 

2 — 71 

\(n ~ 2)a + 2r(a + /3')} - — = 0, (56) 

J 71 + 1 



vr 2-n + r 2( a // + 2j3 „j _ 2r ^ a , + ^ + 2a = 0. (57) 



It is easy to check that the last equation (|57p follows from 
the first 3 equations. This system of equations can be 
reduced to the single third order ODE (master equation) 
for the function v. This can be done by adopting the 
procedure described in the previous section to linearized 
case. Namely, let us find f}' from (I55[) and substitute it 
into (|55|) and (|56|) . Use these two equations to find a 



and a'. Substitute the obtained expression for a' into 
([5]) and differentiate the obtained equation with respect 
to r. After this, again substitute a' in the equation. 
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These operations produce the following master equation 
which contains only v and its derivatives [22| 

nr 2 (y + rv') + 4(n 2 - l)(n + 2)p(v - rv') 
+Ap(n 2 - l)r 2 [rv"' - (n - l)v"] = . (58) 

Before we solve this equation, it proves useful to perform 
a linear change of variables. Let us define the coefficient 
A, with the dimensionality [length] -1 , according to 



V - 



{n 2 - 1) A 2 ' 



(59) 



and introduce a new dimensionless quantity £ = Xr and 
a new dimensionless function, Q(£), according to 



(60) 



Notice that the dimensionality of v is [length] n ~ 2 , while 
Q is dimensionless. The master equation can be cast 
into the form 



+ (4n + 8-£ 2 ) 



; 1)^-^ + 4* + 8)^ 



d£ 2 



0. 



(61) 



The last equation can be rewritten as 



,d 2 T 



dT 



4£ 2 — - 4£(n + 2)— + (4n + 8 - C)T = , (62) 



~^ v " ' 1 di 



where 



Q(0 = 



£(0 



and |=T(0. 



The equation (|6"2")l can be easily solved in terms of the 
modified Bessel functions K q (z) and I q {z). The final so- 
lution for <2(£) has the form 



Q(0 = 



Co 



(63) 



+ £ (n+1)/2 C 1 Kn+ 1 (^/2) + 6*2/^+3(^/2) 



Since we assume that at infinity the spacetime is flat we 
have to put one of the integration constants zero, C 2 = 0. 
In this way we arrive at the final form of the physically 
interesting solution 



Q(0 = ^+& +1)/2 \c 1 K^/2) 
? 2 



(64) 



For completeness, we present also the same solution in 
terms of original variables 



C 



v = — + C 1 A n+1 ^ 2 K n± ,(Xr/2) . 



(65) 



Let us now find the expressions for a and j3. For this 
end we can use the results of the previously considered 
procedure, which involved the equations (JHH) 
and led to the master equation (|58|) . Solving equations 
([55]) . (I56p we find an algebraic solution for a. After 
using (|64j) , we arrive at the following result: 



«(0 



1 



£™n(ri + 1) 
C 



4i 



2 d 2 Q 



d£ 2 

l-n 2 Cl 



(£ 2 + 8)Q 



(66) 



n{n + 1) 

Equation (|55p leads to the equation for (3 



J 2 



d/3 



e (n 
Ci 



1 



2(n + l) 

Integration of this relation gives 



d 2 Q dQ 
d£ 2 d£ 

2K^/2)-^K^/2) 



0(0 



+ 



C 2 



(67) 



^ 2 ±i(^/2)-2n^n^(e/2) 



Substitution of the obtained solutions for «(£)) an d 
/?(£) into the original set of the equations (f55 |) - ((56|) en- 
ables one to make an additional check of the correctness 
of the obtained solutions. 

The modified Bessel functions K^(z), which enter the 
expressions for a and f3, have the following asymptotic 
behavior at the infinity 



1 



4^ 2 



1 



8,-j 



(68) 



Notice, that at the leading order these asymptotics are 
universal, that is they do not depend on the index /x. Us- 
ing (|68p we obtain the following asymptotic expressions 
at large £ for the functions v, a and (3 



X 2 



V, 
C 



Q 



^ + c!^r /2 e-« /2 



n(n + 1) 



1_„ _ gi- n /2 e -S/2 



~ C 2 + 



,72 



(69) 
(70) 
(71) 



VI. FINDING A SOLUTION BY ITERATIONS 

Now we study solutions of the system (I47l) - p9|) with- 
out making an assumption that the gravitation field is 
weak. As we already mentioned in Section 2 the corre- 
sponding black hole solutions has a characteristic scale 
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connected with its gravitational radius tq. For example, 
in the absence of quantum corrections the corresponding 
(Tangherlini) metric is 



ds 2 = -[1 - (ro/r)"- 1 ]^ 2 + 2dvdr + r 2 dw 2 n 
One can always rescale this metric as follows 



ds 2 = r 2 dS 2 



(72) 



(73) 



where dS 2 is the dimcnsionless form of the metric. Mak- 
ing the calculations in this dimensionless form is equiva- 
lent putting ro = 1 in the 'physical' metric ds 2 . In what 
follows we always assume ro = 1. It is convenient to 
rewrite the main system of equations (|47|) - (|49|) in terms 
of the new coordinate x = 1/r. We also introduce the 
new functions M(x) and ty(x) according to 



A = l- x n ~ 1 M(x) , 



* = C. 



(74) 



We denote by a 'dot' a derivative with respect to the coor- 
dinate x. Our normalization means that the asymptotic 
values of M at x — > is always equal to 1. 

In these coordinates and notations the system [jlT])- 
(gSJ) takes the form 



M = 



p(l - n) 



i(x K + xK — K) 



x n - 2 K 2 



-2x n {[{n + l)K + 2xK]M + xKM} 

(75) 
(76) 

K = (n + l)[nxM + 2x 2 M + 2x 2 ^M-2x 3 ^ n ^} 



n 



p{n- l)(n-2) 



x n K 2 



(77) 



In order to solve this system for small p we use the fol- 
lowing procedure. We start with a 'classical' solution. 
Thus we put p — and obtain 



1 , #0 = 0, K = n(n + l)x . 



(78) 



One can replace these expressions into (|75[) and ([76| to 
obtain Mi and ^>i. After this, using (|77|) we obtain K\ 
in which all the terms higher that the first order in p are 
omitted. It is possible to repeat this procedure to obtain 
an approximate solution up to any order in powers of p, 
by means of iterations. 

Denote by {Mk^k, K k } the result of the k-th itera- 
tion. Then one has 



M, 



fc+i 



l+p 



fc+i 



f dxF M {M k ,K k 
Jo 

pF^ k ,K k ), 



(79) 
(80) 



K, 



fe+i 



F { °\M k+1 ^ k+1 )+pF^(K k ). (81) 



Here 



F M (M k ,K k 
(n-1) 



2x n {(n + l)K k M k - 2xK k M k - xK k M k } 



-A(x 2 K k + xK k - K k ) - x n - 2 K 2 k 



[<k] 



2(n-l) _ +1 



[K k - k k v k ] 



[<k] 



Fg\M k ,V k ) = (n + l)[nxM k + 2x 2 M k 
2xH k M k -2x 3 - n ^ k ] 



F { x\K k ) = 



(n - l)(n-2) 



[<k] 



(82) 



In the formulas above the subscript symbol r<k] means 
the expansion into power series up to the order k in the 
parameter p. This prescription provides the correct run 
of the iteration procedure, such that the M k , ^! k and K k 
include powers of p up to the order k. 

The results of the analytical calculations of the first 
two iterations for M, \& and K obtained by using the 
Maple program and checked by using Mathematica are 
given below: 



M, 



1 - [n - 2)(n- l)(n + l)px n+1 
2 



(83) 



(n - 2){n 2 - l) 2 p 2 [{An 2 + 15n + 6)a; 2n ^ 



4(n+l)(2n+ l)x 



n+31 



*2 = 



4p^x 



(n - l) 2 {n + \y{n 2 - 4)(2n + 1), (84) 



= n(n + V)x 

- 2{n + l) 2 (n-2)(n-l){2n+l)x n+2 p (85) 

- -(n-2)(n- l) 2 (n + l) 3 p 2 x 
n 

x [ (8n 3 + 69n 2 + 66n + 16) a; 2 ™ +3 

- 8(n + l)(n + 2)(2n+l)x n+4 ] . (86) 

It is easy to see that the solution presented above include 
zero order terms (which are of course the same as in 
(175)1 ). first order and second order terms, all of them 
have corresponding powers of the expansion parameter 
p. It is not difficult to perform further iterations, but 
starting from the third order the formulas become very 
bulky and hence we do not show them. (The third order 
iteration results for M3, ^3, and K3 can be found in the 
Appendix.) 
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n=3, a=.3e-l 




0.2 0.4 0.6 0.8 1 



x 

FIG. 1: Iterations of the mass function M up to the second 
order. The lines 0, 1, and 2 show the functions Mq(x), M\{x) 
and Mz{x), respectively. The plot is constructed for n = 3 
and a = 0.03. 



n=3, a=.3e-l 




x 



FIG. 2: Iterations of K/Ko up to the second order. Hereto = 
n(n + l)x is the value of K for p = 0. The line is zero order 
in p result. Lines 1 and 2 show the functions Ki(x) / Kq(x) 
and K2(x)/Ko{x), respectively. The plot is constructed for 
n — 3 and a — 0.03. 



n=3, a=.3e-l 



0.04- j 
0.03- 

C 

0.02- 
0.01 1 

o \-:-r--i-rr^--T-rr--rr--r:'' 

0.2 0.4 0.6 OS 1 
x 

FIG. 3: Iterations of C(x) up to the second order. The func- 
tions Co and C\ vanish, so that the plot shows only C2 func- 
tion. 



The plot in Figure Q] shows that the mass function 
M(x) grows from its value 1 at infinity to a larger value 
at the horizon. A natural interpretation of this result is 
that introducing the the Weyl term to the action is qual- 
itatively equivalent to the addition of a negative mass 
density distribution into the black hole exterior. Let us 
remember that we assume that the Weyl term in the ac- 
tion produces the change of the black hole solution which 
is close to the one of the vacuum polarization effect. 

We use equation 

M(x„) = x x H n , (87) 

to obtain a position of the horizon. We also use the 
expression (I43[) for the surface gravity 

k = -- [x 2 e c A] x=XH . (88) 

The results of the calculations of xh and n are given in 
the Table. 



a 




0.01 


0.02 


0.03 


XH 


n — 3 


0.994 


0.987 


0.979 


XH 


n = 4 


0.987 


0.970 


0.950 


k/k 


n = 3 


1.014 


1.028 


1.041 


k/k 


n = 4 


1.036 


1.074 


1.124 



TABLE I: Position of the horizon, xh, and the surface gravity 
k of distorted black holes. 



This table shows that, as expected, for smaller values 
of a, the solution are closer to their unperturbed ones. 
The size of the event horizon, tq/xh, is always greater 
than its unperturbed value ro- As far as the value of a 
increase, the gravitational radius r# = xjj 1 grows up and 
the surface gravity grows up as well. 



VII. 4D CASE 

As we already mentioned in the Introduction, the AD- 
case is a special one. Let us consider it in more details, 
using the results of the iteration procedure presented in 
the previous section. 

The iteration formulas f55|). (fgl)) and (J55J) imply that in 
the 4-dimensional case, when n = D — 2 = 2, the correc- 
tions to the 'classical' solution ([75)1 vanish. One can check 
this directly by substituting (|75|) into the equations (|75|1 - 
(|77p . In principle, there might exist non-perturbative 
in p solutions close to the Schwarzschild metric. Let us 
demonstrate that at least for small perturbations which 
can be treated in the linearized approach, such regular 
both at the infinity and at the horizon solutions do not 



9 



exist. For this purpose let us consider small perturba- 
tions near the 'classical' solution 

M = 1 + em(x) , * = ev(x) , K = &x [l + ek(x)] . (89) 

Here e is a small parameter which should be set to unity 
after performing the series expansion. Substituting these 
relations into the basic equations (jTS"|) - ([??)) one obtains 
that the terms of the order of e° vanish. Keeping leading 
linear in e terms one obtains a following set of equations 



6px 2 [2x(l - x)k + (6 - 7x)k + 3fc] 



(1 - 6px 3 



I8px m = , 



6px 3 [xk + 2k - v] = 0, 



m + xfn — k — (1 — x)v = . 



(90) 



(91) 



(92) 



We obtain now a master equation for this system. For 
this purpose we use (|9ip to find v(x) 



6px 3 



6px 3 



- (xk + 2k) 



(93) 



Let us substitute this expression into (|9"2"j) . Using the 
obtained equation and ([90]) one can solve these linear 
with respect to m and rh equations to obtain the both 
quantities. One gets 



= k 



6px 3 



12px 3 
6px 2 



(1 - 6px 3 )(l + I2px 3 ) 
+ (7x — 6px* — 6)k 



x(l - x)k + (4 - 5x)k , (94) 

: — — — ; — — — 2x(x - l)(ipx 3 + l)k 

(95) 

Differentiating (|94[) and putting it equal to m defined 
by (|95p one obtains a master equation for k. From the 
structure of the expressions for m and m it is evident 
that the master equation does not contain k, but only 
its derivatives. This means that k = const is a solution. 
This solution results in a simple change of the gravita- 
tional radius ro which we originally put equal to 1 . We do 
not consider this renormalization ambiguity and in what 
follows will keep ro = 1. 

Let us now study possible non-trivial perturbed solu- 
tions. By denoting k = Y we write the master equation 
in the form 



F 2 Y + F 1 Y + F Y = 0, 



(96) 



= 

Fi = 
F 2 = 



108p 2 x 6 
- 27px 4 



1 + 108px 2 - lUpx 3 - 
12px 3 (5-6x + 21px 3 

-lUp 2 x 6 + 216pV), 
6px 4 (l - x)(l + I2px 3 ){l - 6px 3 ) 



3„9 



1728p J x 



(97) 



As any second order ODE, can be written in a 
self-adjoint form 



dx (j 1 dx ) ^° 
By comparing and (|9"8|) we have 

fx _ Fi fo _ Fq 
fi F 2 ' /i F 2 

Integrating the first of these equations we obtain 
(l-x)V° 



h = 



1 - 6px 3 )(l + 12pa; 3 



(98) 



(99) 



(100) 



A general solution contains as a common factor an arbi- 
trary integration constant C. We put this constant equal 
to 1. With this choice and p > —1/12, fi(x) is positive 
in the black hole exterior. Using the second equation of 
(IM)) we get 



fo 



{l-x)x 6 Q 



6p{l - Qpx 3 ) 2 (l + 12px 3 ) 2 ' 



(101) 



Q = 1 + I08px z - lUpx 3 - 108p z x b + 1728p J x y .(102) 

The horizon, x = 1, and the infinity, x — 0, are singular 
points of the equation (|9"51) . Let us consider first this 
equation in the vicinity of the horizon. Denote x = l — y, 
then in the region y m one has 



d 2 Y ^dY 
^ dy 2 dy 



cY = , 



1 - 30p - 



6p(l + 12p) 



(103) 



Putting Y ~ y a one finds that a(l + a) = 0. It means 
that one of the two linearly independent solutions is sin- 
gular, Y ~ y 1 , at this point. Hence only one of these 
two solutions is regular at this point, and it is uniquely 
determined by the boundary condition Y{y = 0) = Yq. 

Let us analyze solutions of (|98p near the infinity, that 
is near x = 0. In the close vicinity of this point one has 



F 2 - 6px 4 , F x ~ 60px 3 , F - 1 . 
and the equation (|96p takes the form 

6px 4 Y + 60px 3 Y + Y = . 



(104) 



(105) 



Changing the coordinates r = 1/x and keeping leading 
at r — > oo terms we obtain the following asymptotic form 
of the equation 



6p 



d 2 Y 



dr 2 



— + rY = . 
dr 



(106) 



The asymptotic form of the solutions at infinity is 

Y ~ exp(±r/2A) , A = v/^6p . (107) 
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FIG. 4: Function Q(x,p) in the domain x £ [0,1], p 6 
[0, —1/12]. It is clear from this plot that this function is pos- 
itive everywhere in this domain. 



A solution of (|106[1 which decreases at infinity is 

Y = C*(r 4 + lOAr 3 + 45A 2 r 2 + 105A 3 r + 105A 4 )e~ r/A . 

(108) 

Since a regular solution is uniquely fixed by its value 
at the horizon, in a general case such a solution would 
become increasing at infinity. Thus in a general case 
Y = is the only solution which is regular both, at the 
horizon and infinity. Based on this analysis, one however 
cannot exclude the existence of regular solutions for some 
discrete values of p. We demonstrate now that it never 
happens when p > —1/12. 

For this purpose let us multiply (198|) by Y and integrate 
the obtained expression in the interval x G (0, 1). After 
integration by parts one has 

f dx[hY 2 - f Y 2 } = [hYY]l . (109) 
Jo 



VIII. DISCUSSIONS 

We study black hole solutions in the theory with the 
Weyl action correction to the Einstein gravity in higher 
dimensional spacetimes. We demonstrated that in the 
higher dimensional case with D > 5 this correction re- 
sults in the modification of the Tangherlini metric, which 
is uniquely determined the one parameter, the value of 
the gravitational radius. We developed the iteration pro- 
cedure which for a small value of a = ajr\ allows one to 
obtain a solution outside the horizon. This iteration pro- 
cedure does not work uniformly. In particular for a given 
value of a there always exist such a value of r where the 
iteration is not convergent. For small enough value of a 
this value of r is inside the gravitational radius. In this 
domain terms with higher in curvature corrections must 
play an important role. It is an interesting problem to 
develop a method of solving the equations in this domain 
either by obtaining a reliable analytic approximation or 
by developing stable numerical schemes. Another inter- 
esting open question is to study uniqueness of the exterior 
solution obtained by the iteration method. We analyze 
a similar problem in the 4D-case where the situation is 
much simpler. We first demonstrate that the classical 
Schwarzschild metric is a fixed point of the iteration pro- 
cedure. After this we studied metrics in the vicinity of 
the Schwarzschild metric and demonstrated that at least 
for the value of a < 1/2 (p > —1/12) such regular per- 
turbative solutions are absent. We also showed that in 
a general case additional to the Schwarzschild solutions 
do not exist. It is interesting to generalize these results 
to the higher dimensional case. It should be emphasized, 
that the action with the C 2 term does not modify the 
AD black hole solution. However a similar to the higher 
dimensional case happens, for example, if one considers 
a correction to the Einstein action of the form /(C 2 ), or 
even more general functions of the curvature invariants. 



For a solution which decreases at infinity and is bounded 
at the horizon x = 1, where /i(l) = 0, the expression in 
the right-hand side vanishes. For p < —1/12 the func- 
tion fi(x) is positive in the interval x € (0, 1). Numeri- 
cal analysis shows (see Figure 2]) that for p < —1/12 the 
function fo(x) is negative in the interval (0,1). Hence 
the integrand in the left hand side is non-negative every- 
where in the interval x € (0,1). Hence both Y and Y 
must vanish identically. This proves that for p < —1/12 
the equation has only one bounded solution, Y(x) = 0. 
It means that for small p the Schwarzschild metric is the 
only one solution of the modified by C 2 correction Ein- 
stein gravitational equations (at least in the close vicinity 
of this solution). 
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Appendix. Results of iterations in the third order 

Here we present the third order results for the iteration 
procedure developed in the section 6. 

I 



M 3 = 1- (n-2)(n-l)(n+l)px" +i -- (n-2)(n -1)V {An z +lhn+Q)x Zn+z -A{n+l){2n+l)x 



n+3 



4(n- 2)(n 2 - 1) 
3^ 



(68n 4 + 745n 3 + 1699n 2 + 1166n + 240)x 3 ™ +3 



-12(n + l) 2 (12n 2 + 97n + U)x Zn+4 + 48(n + l)(n + 2)(n + 3)(2n + l)a; 



n+5 



(Al) 



_ 4(n-l) 2 (n + l) 3 (n 2 -4)(2n + l) 
n 



p x 



+ 



32(n-2)(n- l) 3 (n + l) 4 3 



(4n 4 + 41n 3 + 86n 2 + 58n+12)a; 3 ™ +2 - (n + 2)(n + 3)(n + 4)(2n+l)a; 2 " +3 . (A2) 



i4T 3 = n(n + l)x - 2(n + l) 2 (n-2)(n-l)(2n+l)px 



n+2 



2(n-2)(n-l) 2 (n + l) 3 2 



n 



8n 3 + 69n 2 + 66n + 16) x 2n+3 - 8(n + l)(n + 2)(2n + l)x 



n+4 



8(n-2)(n-l) 3 (n + l) 4 3 
3^ P 



48(n + 3) 2 (n + l)(n + 2)(2n + l)x 



n+6 



6(n + l)(32n 4 + 399n 3 + 1005n 2 + 828n + 208)a; 



2n+5 



+ (142n 5 + 1643n 4 + 5222n 3 + 5878n 2 + 2700n + 432)a; 3 " +4 

I 
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